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Abstract 

We discuss some of the properties of the 'collision' of a quantum mechanical 
wave packet with an infinitely high potential barrier, focusing on novel aspects 
such as the detailed time-dependence of the momentum-space probability density 
and the time variation of the uncertainty principle product Axt ■ Apt- We make 
explicit use of Gaussian-like wave packets in the analysis, but also comment on 
other general forms. 



I. Introduction 

The use of wave packets to analyze the non-trivial time-dependence of quantum 
mechanical systems is one important aspect of the study of the classical-quantum 
interface. Popular simulation packages [1] can help students visualize the evolution of 
quantum states (as opposed to time-independent stationary state solutions [2] , [3] , [4] , 
[5]) by allowing them to continuously change parameters (such as the initial width of 
a wave packet) to study what effect they have on the system under study. 

A number of authors have considered various one-dimensional quantum mechanical 
problems in a wave packet approach studying transmission and reflection from square 
barriers [6], [7], [8], [9] or linear potential steps [10], bound state wave packets in sin- 
gle square wells in either position space [11] or momentum space [12], in double wells 
[13], [14], or in systems of relevance to solid state physics [15], [16]. Such examples 
of wave packet behavior are also increasingly useful as teaching tools since the behav- 
ior of Coulomb wave packets on circular [17] or elliptical [18] orbits are being tested 
experimentally on Rydberg atom systems [19], [20]. 

Numerical methods for solving the time-dependent Schrodinger equation have been 
discussed [6], [21] and, in special cases, closed-form analytic results can be obtained by 
use of the time-development operator [22], [23], e~ lHt ^ h ^p(x,0) = ip(x,t), to solve the 
initial value problem. Another approach is to combine a large number of individual 
stationary state solutions for both unbound [24] and bound state problems [11] to 
obtain wave packets. The most familiar example is the explicit calculation of the 
Gaussian free-particle wave packet which is treated in the majority of elementary texts. 
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A simple variation is to consider a particle, subject to the one-dimensional 'infinite wall' 
potential [5] given by 



so that it is free for x < 0, but a wavepacket impinging on the 'wall' at the origin 
will 'bounce'. Andrews [25] has shown how some of the most obvious aspects of the 
'collision' process (namely the long-time development of the reflected wave packet, 
interference effects during the 'collision', etc.) can be understood by considering com- 
binations of free-particle 'mirror' solutions, and we will use some of his arguments. 

In this note we will examine, in some detail, the 'bounce' of free-particle wave 
packets from the infinite wall potential described by Eqn. (1), focusing on several 
other issues, namely the behavior of the momentum-space wave packet solutions, the 
widths of the position- and momentum-space packets during the 'bounce', and the 
uncertainty principle product Ax t ■ Ap t as a function of time. We will make extensive 
use of the free-particle Gaussian wave packet in our discussion, but we also present 
results for other, more general forms; for completeness sake, however, we very briefly 
review the essentials of the Gaussian case. 

A free-particle wave packet can be constructed, using any initial momentum-space 
weighting function, 0), via 



to give a time-dependent position-space wavefunction, ip(x,t). The momentum-space 
solution itself has a trivial time-dependence, namely 




for x < 
for x > 



(1) 





4>[p,t) = 4>[p, 0)e 



ip 2 t/2mh 



(3) 
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For the case of a Gaussian momentum-distribution, here written in the form 



0i (p, 0) = l® e -« 2 (P-Po) 2 /2 e ~ipxo/n (4) 
V 

the necessary integral in Eqn. (2) can be done to obtain the well-known result 

ijj^X = - c i\po(x-xo)-p2t/2m\/h c -(x-x ~p t/m) 2 /2a 2 K 2 F ^ 

\JahFyfjr 

where F — l + it/t and t = mha 2 . This solution describes a Gaussian position-space 
wave packet whose width increases with time, characterized by arbitrary initial values 
of xq and po- The resulting position-space probability density is 

Pfreefrt) = |^l(x,t)| 2 = -1= e"**"* o-PO*M) a /^ ( 6 ) 

where (3 t = + (t/to) 2 ] 1 ^ 2 and various important expectation values are given by 



(x) t = x + ^- , Ax t = ^Jl + (Pl , (p) t =p , &Pt = (7) 



m 2 V Vto/ v2a 

For the infinite wall case, we also can obtain wave packet solutions from Eqn. (2) 
by substituting the appropriate plane wave solutions 



ipx/n , e **/ft- e -**/ft for:r<0 

^ lor ,r > 



in the basic integral. In this approach, the integrals must be performed numerically. 
On the other hand, we can also make use of the method of Andrews [25] and use any 
free-particle wave packet solution ip(x,t) via 

which satisfies the Schrodinger equation for the potential in Eqn. (1) as well as the 
appropriate boundary condition at the wall. In either case, if the original free-particle 



wave packet is properly normalized, the 'bouncing' wavepackets will also be very close 
to being normalized, provided they are initially far enough from the wall so that any 
contribution from the 'tail' in the x > region is negligible. In either case, how- 
ever, in order to obtain the time- dependent momentum-space wavefunction, we must 
numerically evaluate the Fourier transform 

<f>(p, t) = - J= £J e -*»/» i/;(x, t) dx (10) 

To illustrate the behavior of such a 'bouncing' wave packet, we show in Fig. 1 plots 
of the position- and momentum-space probability densities for a Gaussian wavepacket 
for various times before and after a collision. We have used the following values in 
numerical integrals: 

H — 1 , m — 1 , po = 10 , x = —10 , a — 1 (11) 

With these values, the spreading time is t — 1 and the time it takes the packet to 
return to its initial starting point is T = 2t = 2, so that an appreciable amount of 
spreading is obvious. In order to see what features of such 'collisions' are specific to 
Gaussian packets, in Fig. 2 we show the same plots, but for an initial momentum-space 
amplitude given by a Lorentzian form, namely 

MP, 0) = J^-T^f ~ vTXTT e ~ WX ° /h ( 12 ) 

V 7T [a 2 (p -p ) + !] 

The corresponding initial position-space wavefunction is 

^2(^,0) = _L e -l*-*o|/fi» ( 13 ) 
V ah 



but the further time-dependence can only be evaluated numerically using Eqn. (2). 
Using these two cases, we can make some general comments: 



(i) The non-Gaussian position-space wave packet comes to approach the Gaussian 
form more and more closely, as it evolves in time. This behavior is seen for a 
large number of other, single- humped initial distributions [26]. 

(ii) The momentum-space probability density well after the collision is related to 
the initial density by P a f ter (p,t) = Pbefore(—p,t) corresponding to the reversal of 
each momentum component during to the collision. 

(hi) At the moment of the collision, however, the momentum distribution is not sym- 
metric. This is clearly due to the fact that the high momentum components are 
preferentially in the leading edge, and are the first to be reflected to negative 
values. This also implies, as will be seen later, that the expectation value (p) t 
is slightly negative at t — Tc, the collision time. We note that other 'velocity 
effects' have been discussed [7], [8] for various kinds of wave packet scattering. 

To focus on the details of the collision event, in Fig. 3 we plot the time-dependent 
\ip(x, t)\ 2 and \(j)(p, t)\ 2 (for the Gaussian wave packet) for times nearer the actual 
'bounce', bracketing t = T c , and we note some additional aspects of the process: 

(i) The time-dependence of <p(p,t), which is non-trivial only during the collision, is 
more clearly visible as is the eventual return to 'symmetry' of \4>(p, t)\ 2 . 

(ii) The spread in the position-space probability density at the time of the collision 
is substantially smaller than Ax t either immediately before or after the collision. 
We will address this point below, using an analytical evaluation of Ax t =T c - 

In order to examine more of the differences between the purely classical and quan- 
tum approaches to the collision of a point particle, we plot in Fig. 4 calculations of 

5 



the expectation values (x) t and (p) t , which are easily evaluated numerically. In this 
figure, we show the expectation values of x and p (solid curves) for the bouncing wave 
packet, as well as those for the free-particle wave packet (dotted curves) for the stan- 
dard set of parameters in Eqn. (11) except that we use the value of a — 0.5; this value 
is chosen to make the spreading of the wave packets more obvious since in this case 
t = mha 2 is much smaller. We also indicate the 'one sigma' limits given by (x) t ±Ax t 
and (jp) t ± Ap t as dashed curves. 

We first note that the guaranteed relationship between (x) t and {p) t , namely {p) t = 
md{x)t/dt, is trivially observed long before and long after the collision, while the same 
qualitative connection between (x) t and (p) t near t = Tq is now also apparent and 
different than a purely classical 'bounce' for a point particle which would have a cusp 
(discontinuity) in x(t) (v(t)) at the collision time. Finally, we can see that the position 
spread at the collision time, Ax t =T c is slightly smaller for the 'bouncing' wavepacket 
than for the free-particle packet with the same initial parameters, while the momentum 
spread is much larger at T c than in the free-particle case. 

In order to examine the 'compression' of ip(x,T c ) and related issues, we plot in 
Fig. 5 the values of Ax t , Ap t and the uncertainty principle product Ax t • Ap t (in units 
of h) for a range of values of a, but keeping the other parameters fixed as in Eqn. (11). 
In each case we see that Ax t=Tc is indeed smaller than its value for the free-particle 
wave packet and in the cases where a > 1, it is even smaller than it's original spread, 
Ax t =o. This effect is perhaps intuitively obvious as the high momentum components 
are reflected first, while the low momentum pieces 'pile up', leaving the position-space 
wave packet temporarily narrower. This is not a violation of the x — p uncertainty 
principle as many other cases of such behavior are known; for example, similar effects 
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are seen in explicit constructions of wave packet solutions [27] or more simply in the 
direct examination of the time-dependence of the uncertainties in x and p in complete 
generality [28], for the harmonic oscillator potential. In our case, the fact that Ap t 
does indeed increase during the collision is even more obvious, especially from the 
time-dependence of \(f>(p, t)\ 2 shown in Figs. 1 and 2. During the collision, instead of 
being dominated by the intrinsic width of a single (f>(p, t) peak, Ap is dominated by 
the distance between the peaks. 

^From the explicit numerical calculations used to generate Fig. 5(a), we find to an 
excellent approximation that the position-space spread at the collision time is given by 



and we can make use of the more analytic approach followed by Andrews, at least at t — 
Tc where the expressions simplify dramatically, to understand this effect quantitatively. 
Using the explicit ipi(x, t) in Eqn. (5) and the expression in Eqn. (9), we can construct 
an excellent approximation to the 'bouncing' wave packet for the Gaussian case, namely 



which is approximately normalized, and the error is exponentially small for the pa- 
rameters we use, namely po!3t c /Ti >> 1- For these values, the sin 2 (p Q x/h) variation 
can very reasonably be replaced by its average value of 1/2 and the resulting integrals 
performed exactly. We then have, to an excellent approximation, 




0.60 



(14) 




(15) 



(x)t c 




and 




(16) 



so that 



all of which are observed numerically! 

A similar semi-analytic result can be obtained which describes the expectation 
value of the momentum operator at the collision time, namely {p)r c - The values of 
(p)t required for Figs. 4 and 5 have been obtained numerically, by using the momentum 
space probability density, but using the wavefunction representation in Eqn. (9), we 
can also obtain an explicit formula for (p) t for the special case of t — Tc using the 
representation of p as a differential operator acting on the position-space wavefunction. 
For the free-particle wave packet we naturally have 



/+oo 
^i(a;,f)*p^i(a;,t) dx = 
-oo 



Po 



(18) 



and the evaluation is straightforward and well-defined for all times. In contrast to this 
case, if we naively attempt to evaluate {ip\p\ip) in this way we find that the expectation 
values are not necessarily Hermitian due to the 'asymmetry' caused by the presence 
of the wall. If, however, we instead adopt the 'symmetrized' version (which reduces to 
the standard value for the free-particle case) 

1 r° 



t c = 2 J [(p^( x , T c)) i/>(x, T c ) + ip*{x, T c ) (pip{x, T c 



dx 



(19) 



we find that (p) is guaranteed to be real. Using this trick, we find that the expectation 
value at the collision time is 



Ah t 



7rW A 3 



Tc/to 



[ ,c siir ( ) < '' ■ dx 

J — oc 



. 2 (PqX\ - x 2/& 
sin I I p 1 

V h 



(20) 



/1 + (Tc/to) 2 _ 

where sin 2 (p ^/^) term is replaced by its average value of 1/2. This analytic approxi- 
mation agrees with all of our numerical calculations to the desired accuracy. 
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Figure Captions 



Fig. 1. Plots of a Gaussian wave packet striking an infinite wall (bold vertical line.) 
\ip(x,t)\ 2 versus x is shown on the left, while the corresponding \(f>(p, t)\ 2 versus 
p plot is shown on the right. The initial position x is shown on the left, while 
the values of the central momentum long before (+Po) and long after (— p ) the 
collision are indicated on the right. The numerical values used are those in 
Eqn. (11). 

Fig. 2. Same as Fig. 1, but for a wave packet described by an initial Lorentzian (given 
by Eqns. (12) and (13) ) which then evolve in time. 

Fig. 3. Same as Fig. 1, but for times nearer the actual 'collision' at t ps T c . Note that 
the momentum distribution at the moment of collision is not symmetric. The 
peak near +p Q is skewed towards values with p < +po, while the feature near 
— p is similarly enhanced with values just below — p . 

Fig. 4. Plot of (x) t (top) and (p) t (bottom) as a function of time over the same time 
interval as shown in Fig. 1. In order to emphasize the spreading of the wave 
packet, we use a = 0.5; otherwise the parameters are as in Eqn. (11). Results 
are shown for the bouncing (solid curves) and the free-particle (dotted curves) 
packets. The dashed lines indicate the one standard deviation spreads in each 
case, calculated numerically. Note that the expectation value of momentum at 
the collision time t = T c is slightly negative. 

Fig. 5. Plots of (a) Ax t , (b) Ap t and (c) the uncertainty principle product Ax t ■ Ap t 
(in units of Ti) versus t for various 'bouncing' wavepackets over the same time 
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interval as in Fig. 1. The various cases corresponding to a = 3, 2, 1, 1/2, 1/3 are 
given by the dash- dash- dot, dash, solid, dot-dash, and dotted curves. Otherwise, 
the standard set of parameters in Eqn. (11) are used. All the Gaussian packets 
shown start with Ax t ■ Ap t = h/2. 
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